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In many integrable models static (equal time) correlation functions of local observables after a 
quantum quench relax to stationary values, which are described by a generalized Gibbs ensemble 
(GGE). Here we establish that the same holds true for dynamic (non equal time) correlation func- 
tions. More generally we show that in the absence of long-range interactions in the final Hamiltonian, 
the dynamics is determined by the same ensemble that describes static correlations. When the latter 
is a GGE the basic form of the fluctuation dissipation theorem holds, although the absorption and 
emission spectra are not simply related as in the thermal case. For quenches in the transverse field 
Ising chain (TFIC) we derive explicit expressions for the time evolution of dynamic order parameter 
correlators after a quench. 



Introduction. By virtue of their weak coupling to the 
environment ultra-cold atomic gases provide ideal test- 
ing grounds for studying nonequilibrium dynamics in iso- 
lated many-particle quantum systems. A number of re- 
cent experiments [TH5] have observed essentially unitary 
time evolution on long time scales. This has stimulated 
much theoretical research on fundamental questions such 
as whether observables generically relax to time indepen- 
dent values, and if they do, what principles determine 
their stationary properties. Relaxational behaviour at 
first may appear surprising, because unitary time evolu- 
tion maintains the system in a pure state at all times. 
However, it can be understood intuitively as a property 
of a given finite subsystem in the thermodynamic limit, 
with the role of the bath being played by the rest of the 
system. 

Dimensionality and conservation laws strongly affect 
the out-of-equilibrium dynamics. Ground breaking ex- 
periments by Kinoshita, Wenger and Weiss [5] on trapped 
87 Rb atoms established that three dimensional conden- 
sates "thermalize" rapidly, i.e. relax quickly to a sta- 
tionary state characterized by an effective temperature, 
whereas the relaxation of quasi one-dimensional systems 
is slow and towards an unusual non-thermal distribution. 
This difference has been attributed to the presence of 
additional (approximate) conservation laws in the quasi- 
1D case, which are argued to constrain the dynamics 
in analogy with classical integrable systems. The find- 
ings of Ref. [5] sparked a tremendous theoretical effort 
aimed at clarifying the effects of quantum integrability 
on the non-equilibrium evolution in many-particle quan- 
tum systems, see e.g. Refs [7H22] and references therein. 
A widely held belief that has emerged from these stud- 
ies is that the reduced density matrix of any subsystem 
(which determines correlation functions of all local ob- 
servables within the subsystem) is described in terms of 
either an effective thermal (Gibbs) distribution or a so- 
called generalized Gibbs ensemble (GGE) [8]. The former 
is believed to represent the generic case, while substantial 
evidence has accumulated that indicates that the latter 
arises for integrable models. 

Given this state of affairs, a natural question to ask is 



whether the GGE applies only to static properties in the 
stationary state, or whether it holds more generally even 
for dynamical properties at late times after the quench. 
The latter would open up the possibility to explore the 
physics of the GGE in its full glory through experimen- 
tal probes at finite energies, such as photoemission spec- 
troscopy [23]. In the first part of this letter we prove 
quite generally, that dynamical correlations of local op- 
erators acting within a given subsystem in the stationary 
state after a quantum quench are determined by the same 
distribution function as static correlations. In particular 
this means that whenever the GGE describes static cor- 
relations in the stationary state, it also applies to the 
dynamics. 

Stationary State Dynamics after a Quantum Quench. 
We consider the following quench protocol. The system is 
prepared in the ground state |\&o) of a Hamiltonian H(h n ) 
with local interactions, where h is a system parameter 
such as a magnetic field. At time t = we suddenly 
change ho to h and the system time evolves unitarily 
with Hamiltonian H(h) thereafter. We are interested in 
expectation values of the form (ti, . . . , t n > 0) 

(* (*)|Ol(*l)...On(*»)l*o(t)), (1) 

where Oj are local observables. We wish to demonstrate 
the following. If the stationary state of a quantum many- 
body system after a quantum quench is described by a 
density matrix p s tat such that for observables Oj acting 
only within a subsystem S one has 

lim (Mt)\Oi . . • On|*o(*)> = ^(fttatOi ...O n ), (2) 

then dynamical correlations are described by the same 
density matrix, i.e. for t±, . . . , t n fixed we have 

lim (* (*)|Ol(*l)...On(*n)l*0(*)> 
t—too 

= Tr(p stat 1 (t 1 )...O n (t n )) . (3) 

The proof of this statement is based on the Lieb- 
Robinson bound [24] and more specifically the following 
theorem by Bravyi, Hastings and Verstraete [35] : let Oa 



2 



be an operator that differs from the identity only within 
a local region A. Now define the projection of the (non- 
local) operator Oa(*) to the subsystem S D A by 



(S ) _ trg[Q A ft)]®Ig 
( } = trg[I g ] ' 



(4) 



where S is the complement of S. If the time evolution is 
induced by a short-range Hamiltonian, then 

\\0 A {t)-Of (f)||<c|A| e -^, (5) 

where ||.|| is the operator norm, v is the maximal velocity 
at which information propagates [24 , d is the (smallest) 
distance between S and A, \ A\ is the number of vertices in 
set A, and £, c positive constants. Assuming the operator 
02 to be bounded, ||02|| < we therefore have 

K*0i(ti)0 2 (t2)}t| < ||«M*i)0 2 (t 2 )|| 

< \\SOi(h)\\ \\0 2 (t 2 )\\ < Ci \Ax\k e -^^ , (6 ) 



where 



denotes expectation value with respect to 



|*o W> and 5£>i(t) = O x {t) - Of(t). The first inequality 
holds because the operator norm is an upper bound for 
the expectation value on any state, while in the last step 
we used Eqn ^ implies that 



([[O^t^t = (Of(ti)0 2 (t2))t + ai(ti,t 2 )e— "V" , 

(7) 

where ai(ii,i 2 ) is a bounded function. By repeat- 
ing the steps leading to §f§ for the operators 2 {t 2 ) 
we arrive at (0 1 (t 1 )0 2 {t 2 )) t = (Of (ti)0f(t 2 )) t + 
Si=i a i(^ii ^2) exp (— ) ; where a 2 (ii, t 2 ) is another 

bounded function. We may now use the assumption |2]) 
for the expectation value on the right hand side since all 
operators act within subsystem S 



lim (0 1 (t 1 )0 2 (t 2 )) t = Tr(p stat OS(h)Ol(t 2 ))- 



a i (ti 1 t 2 )e « 

i=l 



(8) 



The chain of inequalities ([6]) also holds for the average 
with respect to the density matrix p s tat> i-e. 



Tr(p Btat 1 (t 1 )0 2 (t 2 )) = Tr(p stat Of (f x )Of (t 2 ))- 
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E 



bi(tx,t 2 )e 



(9) 



where bi(t\,t 2 ) are bounded functions of ti, 2 . Finally, 
combining ([8| and Q and then taking the size of the 
subsystem S to be infinite we obtain Q in the case n = 2. 
The generalization to arbitrary n is straighforward. 



Generalized Gibbs ensemble. We now concentrate on 
a quantum quench in an integrable model in one dimen- 
sion with Hamiltonian H (h) and local conservation laws 
I n , i.e. [I m ,I n ] = [I mi H(h)\ = 0. The full (reduced) 
density matrix of the system (of a subsystem A) at time 
t after the quench is 

p(t) = |* (t)X*o(*)|, pa(*) = 1^(p(*)). (10) 

where A is the complement of A. It is widely believed, 
and was shown for quenches of the transverse field in the 
TFIC in Refs HE 123, that 



lim p A (t) = Tr^(pGGE) 



where 



1 



Pgge — 



(11) 



(12) 



■^GGE 

is the density matrix of the GGE and _Zqge ensures 
the normalization tr(pGGE) = 1- Eqn ( |ll| establishes 
that all local, equal time correlation functions of a given 
subsystem in the stationary state are determined by the 
GGE (12 1. Applying our result ^ to the case at hand, 
we conclude that dynamic correlation functions are also 
given by the GGE, i.e. 



lim (#o(*)|Oi(ti)...0n(tn)l*o(t)> 

t— >oo 

= Tr(p GGE 1 (t 1 )...O n (t n )) ■ 



(13) 



Fluctuation Dissipation Theorem. A key question re- 
garding dynamical properties in the stationary state after 
a quantum quench is whether a fluctuation dissipation re- 
lation holds [28]. Given the result ( 13 ) we can answer this 



question for cases where the stationary state is described 
by a GGE. The linear response function of observables 
Aj and Bi acting on sites j and I of a translationally 
invariant lattice of L sites is 



XAs(w,q) = - 



3,1 



dre 1 



-tq(r^-ri) 



Sab(u,ci) 



L ^ 

3,1 



(15) 



x tr[p GGE [A,-(T),B z ]]. (14) 

On the other hand, the spectral function of the same two 
observables in the stationary state is given by 

-1 />oo J 

Loo 2^ 

X tt[pGGEAi(T)Bj] . 

Using a Lehmann representation in terms of Hamiltonian 
eigenstates it is straightforward to show that 

- -Im xab{u, q) = S A b(u, q) - S B a(—u, -q) > (16) 

7T 

i.e. the basic form of the FDT holds. However, as 
was already noted in Ref. [28] for the TFIC, unlike in 
the thermal (Gibbs) case, the negative frequency part 
Sba{—^, — q) is not related to the positive frequency 
part by a simple relation of the form Sab(— — q) = 
q), where /(w) is independent of A and B. 
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Transverse Field Ising Chain. We now focus on the 
dynamics after a quantum quench in a particular exam- 
ple, the TFIC described by the Hamiltonian 



H(h) 



ha* 



(17) 



where er" are the Pauli matrices at site j, J > and 
we impose periodic boundary conditions cr£ +1 = oy\ 



The model (17 1 is a crucial paradigm of quantum crit- 
ical behaviour and quantum phase transitions |29j . At 
zero temperature and in the thermodynamic limit it ex- 
hibits ferromagnetic (h < 1) and paramagnetic (h > 1) 
phases, separated by a quantum critical point at h c = 1. 
For h < 1 and L — > oo there arc two degenerate ground 
states. Spontaneous symmetry breaking selects a unique 
ground state, in which spins align along the ^-direction. 
On the other hand, for magnetic fields h > 1 the ground 
state is non-degenerate and as the magnetic field h is in- 
creased spins align more and more along the z-direction. 
The order parameter for the quantum phase transition is 
the ground state expectation value (cr x ). We note that 



the model (17 1 is (approximately) realized in systems of 
cold Rb atoms confined in an optical lattice [30] . 

Two point dynamical correlation functions are of par- 
ticular importance due to their relationships to response 
functions measured in photoemission and scattering ex- 
periments. The two-point function of transverse spins 
<*o(*)|oj+i(n)o-|(75)|*o(*)> in the TFIC can be calcu- 
lated by elementary means [3T] as it is local in terms of 
Jordan- Wigner fermions. Our goal is to determine the 
dynamical order-parameter two-point function 

after quenching the transverse field at time t = from ho 
to h for times r^a > 0. This can be achieved by employ- 
ing a generalization of the form factor methods recently 
developed in Ref. [35] to the non-equal-time case, and 
augmenting the results obtained in this way by exploit- 
ing the knowledge of exact limiting behaviours derived 
in Refs [26l [32] ■ Our approach is outlined in the Supple- 
mentary Material. For quenches within the ordered phase 
(ho, h < 1) we obtain for large positive £, t 



p xx {£ 7 t + T,i)~ C FF (h , h) R(£, r, t) 



(19) 



where 



R(e, T ,t) 



exp 



77 dk 



log ( cos A k 



C§ F (h ,h) = 



x min { max{4(fc)T, £},e' h (k)(2t + r)} 
1 - hho + y/(l-h 2 )(l-hl) 



2y/l-hhoi/T 



hi 



(20) 



Here £/,(&) = 2JVl + h 2 — 2ft, cos k is the dispersion rela- 
tion of elementary excitations of the Hamiltonian H{K), 



h„=1/3, h=2/3, i=2Q, T/L=16/3 
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FIG. 1: Non-equal-time two point function after a quench in 
the ordered phase from ho = 1/3 to h = 2/3. The distance 
and time T are fixed at I = 20 and T/tp — 16/3 respectively. 

cos Afc = 4J 2 (1 + hh — (h + hO) cos k) / eh(k)eh (k) and 
£ 'h(k) = dsh(k) j 'dk. An important scale in the problem 
is given by the "Fermi-time" 

I 



t F = 



2v„ 



max fc e' h (k) , 



(21) 



where w max is the maximal propagation velocity of the 
elementary excitations of the post-quench Hamiltonian 
H(h). We note that the dominant contribution at large 
£,t (20) has a vanishing imaginary part. This is similar 



to the corresponding correlator at finite temperature in 
equilibrium [29) . In order to assess the accuracy of the 
asymptotic result ( 20 ) at short and intermediate times 



and distances we have computed the correlator ( 19 1 nu- 



merically on large, open chains by means of a determi- 
nant representation and then extrapolated the results to 
the thermodynamic limit. A comparison between ( |19[ ) 
and the numerical results for a quench from h — 1/3 to 
h = 2/3 and distance £ = 20 is shown in Fig. [I] The 
agreement is clearly excellent. The qualitative behaviour 
of p xx (£, T + t,T — t) is as follows: r = corresponds to 
the known [32] equal-time correlator at time T after the 
quench. The correlator remains essentially unchanged 
until t = tp (corresponding to t\ — r 2 = 2tp in (19)), 
where a horizon effect occurs. At later times t > t F the 
correlator decays exponentially. 

For quenches within the disordered phase (ho,h > 1), 
we obtain for i> max (2i + r) > £ 

p xx {£, t+T, t) ~ hCM^o \ h- r )F{£, r, t)R{£, r, t) , (22) 



where R(£,T,t) and Cf F are given by (20) and 



F(£,r,t) = 



* dkJe ilk 
-7T Ke h (k) 



+ 2itan(^)cos(e fe (2t + r))sgn(£-efcT) . (23) 
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h„=2, h=3, £=30, T/t c =16/3 




0.2 



— 0.1 



+ 

I- 

V o 



-0.1 



-0.2, 



h =2, h=3, £=30, T/t F =16/3 



Numerics 
Form Factors 



0.5 




1.5 
rlt 



2.5 



FIG. 2: Real part of the non-equal-time two point function 
after a quench in the disordered phase from ho = 2 to h — 3. 
The distance and time T are fixed at £ — 30 and T/tp = 16/3 
respectively. Data points are numerical results (see the text 
for details) and the solid line is eqn (22 1. 



FIG. 3: Imaginary part of the non-equal-time two point func- 
tion after a quench in the disordered phase from ho — 2 to 
h = 3. The distance and time T are fixed at i = 30 and 
T/tp ~ 16/3 respectively. Data points are numerical results 



(see the text for details) and the solid line is eqn ( 22 1 



In the complementary regime v max (2t + T) < £ the corre- 



lator is exponentially small and the expressions (22 23) 
no longer apply. Outside the "light-cone" v ma x T < ^ the 



first contribution in (23) is exponentially small, whereas 



the second one decays as a power-law. The result (22) 
is obtained by a generalization of the form factor [33] 
approach developed in Ref. (32] and is based on an ex- 
pansion in the density of excitations of H(h) in the ini- 
tial state after the quench. Hence it is most accurate for 
quenches where this density is low and breaks down for 
quenches from/to the quantum critical point. In Figs [2] 
and [3] we compare the asymptotic result ( 19 ) to numer- 
ics obtained in the way described above. The agreement 
for the chosen set of parameters {I = 30, ho = 2, h = 3 
and T/tp = 16/3) is seen to be excellent. The value of 
pxx (£^t + t,T — t) at t = equals the known equal-time 
correlator at time T after the quench [32 , which happens 
to be quite small in the case considered. The correlator 
remains essentially unchanged up to a horizon at r = tp 
(corresponding to t = tp/2 in (191), and for times r > tp 
exhibits an oscillatory t~ x I 2 power- law decay. We note 
that the result (19j 20) for quenches within the ordered 
phase can be obtained in an alternative way by gener- 
alizing the semiclassical approach of Ref. [31] (see also 
Refs [131 [22]) to the non-equal time case and then el- 
evating it using exact limiting results of Refs [3H1 |3"2"] . 
However this method fails to reproduce the result for 
quenches in the disordered phase. 

Having established in the first part of this work that 
the t — > 00 limit of p xx (£,t + Ti,t + T 2 ) is described by the 
GGE, an important question is how quickly this limiting 
behaviour is approached. It follows from |l9 22) that 
for quenches within the ordered (disordered) phase the 



limiting value for fixed t\ % % and I is approached as a t~ 3 
(i~ 3 / 2 ) power law. 

Conclusions. We have considered dynamical (non- 
equal-time) correlation functions of local observables af- 
ter a quantum quench. We have shown quite generally, 
that dynamical correlators of local observables in the sta- 
tionary state are governed by the same ensemble that de- 
scribes static correlations. For quenches in the TFIC this 
implies that they are given by a generalized Gibbs ensem- 
ble, for which the basic form of the fluctuation dissipation 
theorem holds. We have obtained explicit expressions for 
the time evolution of dynamic order parameter correla- 
tors after a quench in the TFIC. 
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Polo fellowship of the University of Bologna and a INFN 
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Supplementary Material 



The dynamical two-point functions (19 22 ) are deter- 



mined by a generalization of the form factor approach 
recently developed in Ref. [55] • The latter is based on a 
Lehmann representation of two-point functions in terms 
of simultaneous eigenstates of the momentum operator 
P and the post-quench Hamiltonian H(h) 

n 

H(h)\kl,...,k n )i = ^2^h(kj) \ki, . . . , fc n ) a , 

J=l 

n 

P\kl,...,kn)i = [^%]|fc 1 ,...,fc„) a , (24) 
J = l 

where a = R, NS correspond to periodic/antiperiodic 
boundary conditions on the Jordan- Wigner fermions [3 2) . 
For a quench within the disordered phase the state 
\^o(t)) has the following representation in a large, but 
finite volume L 



l*o(«)) 



\B(t)) 



NS 



y/ NS (B(t)\B(t)) m 



(25) 



where 



\B(t)) 



NS 



E 



n=0 0<pi,...,p n £NS j=l 

X \-pi,Pl,...,-p n ,Pn)NS, (26) 



K{k) 



sin(fc) (ho — h) 



£h (fc)gfe(fe) 



l^'jf ' + l + hh -(h + h ) cos(fc) 



(27) 



The function K(p) is related to the quantity cos(A p ) de- 
fined in the main text by K(p) = tan(A p /2). The dy- 
namical order parameter two-point function 



p xx (£,t + n ,t + T 2 ) 



NS 



:b\b) 



NS 



(28) 



has the following Lehmann representation 
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NS 



E 



[B(t)\af +m ( n )a* m (T 2 )\B(t)) NS = £ 

J^ s e i(,T 2 -r 1 )£ h (q r )+iq r e 



n\m\ 

m,n=0 °<pi PnENS 

CXfci fc m eNS 



jr(p i )e- ai <*- Ha > eh <«> 



3=1 



2i(t+Tl)E h (fe() 



U=l 



*E E 

s=0 gi,...,g,eR 



ns(5|B) ns =exp[ J! log(l + ^ 2 (?))] • (30) 
o< 9 eNS 



The form factors 



(k nl , k n 



(31) 



are known exactly [33], see eqns (109)-(111) of Ref. 



The leading behaviour of ( 29 1 is evaluated by considering 



it as a formal expansion in powers of the function K{p). 
As shown in Ref. |32j this corresponds to an expansion, 
where the small parameter is the density of excitations 
of the post-quench Hamiltonian H(h) in the initial state 
|$o(0))- We determine the dominant contributions at 
large £, t and \t\ — t%\ to ( p9| for a given order in the 
formal expansion in powers of K(p), and then sum these 
to all orders. The structure of this calculation is similar 
to the equal time case (t% = r 2 ) considered in Ref. [3"2"] . 
but the details differ substantially and will be reported 



(km, -km, ...,kx- fci|0£Jt?i, • • ■ ,q s ) 

x (q s ,...,q 1 \a m \ - pi,pi, . . . , -p n ,p n ) , (29) 

I 

elsewhere. The result of the form factor calculation for 
quenches within the disordered phase is 

p xx (£,t + T,t) ~ 2jVh(h 2 -l)iF{£,T,t)R (l,T,t) , (32) 

where the function F is given above in ( |23[ ) and 

r f* c\k 

R (£,T,t) = cxp -2 / — K 2 (k) 
1 Jo * 

x min|max{4(fc)r,£},e' h (fc)(2t + r)| . (33) 

We now use that the general structure of the resumma- 
tion for ^ is the same as for t\<x = 0. This allows 
us to go beyond the low-density expansion by exploit- 
ing results obtained in Refs [26j G2] for t\ = t% = by 
means of determinant techniques. In this way we arrive 



at eqn (22 1. Quenches within the ordered phase are an- 



alyzed in the same way. 



